Understanding BCS Theory 
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New calculation reveals that E is constant in a thin layer across the Fermi surface, befitting the definition of 
energy gap parameter, A varies dramatically. The BCS self-consistent equation has a simple and exact solution, 
showing that the well-known ratio 3.5 must be replaced by 4. A simple formula is found to estimate energy gap 
from first principles with reasonable accuracy, useful to the current research into high T c superconductivity. 
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According to Bardeen, Schrieffer and Cooper (BCS) 
the superconducting energy gap is measured by the min- 
imum value of E. In the original approximation A is a 
constant that equals the minimum value of E Jjj]. In the 
literature 2A has become synonymous with energy gap. 
However, we find that E is constant in a thin layer across 
the Fermi surface, befitting the definition of energy gap 
parameter, A varies dramatically. Since A has played 
such a pivoting role in the current study of supercon- 
ductivity, this matter appears to be serious enough to 
warrant attention. 

While the BCS theory is deep and fundamental, spe- 
cific results from this theory could be improved, because 
these results arise from approximations employed to ease 
mathematical difficulties. In particular the BCS self- 
consistent equation at T > has a simple and exact solu- 
tion which can be proved by simple inspection. This ex- 
act solution leads to a rather surprising conclusion that 
the well-known BCS ratio 3.5 must be replaced by 4. 

Furthermore, we believe that the BCS theory has sig- 
nificant power of prediction which was unnoticed pre- 
viously. We use the well-established method of itera- 
tion 0] to solve the BCS self-consistent equation. We 
use the Mott- Jones assumption as our first approxima- 
tion H which leads to a simple formula allowing us to 
estimate energy gap quantitatively from first principles, 
the first time to our knowledge. We find from 14 su- 
perconducting metals that on average the energy gap is 
15.7 x 10~ 4 eV, compared with the experimental value 
11.3 x 10~ 4 eV. This could be interesting, in the light of 
the current intensive effort to search how high T c might 
be for various novel superconductive materials. 

We start with the BCS self-consistent equation at 
T > 0: 



E(k + q) 

2k B T 



(1) 



where E = (A + e ) > , e is the electron energy relative 
to the Fermi level, k and q electron and phonon vec- 
tors, respectively jlj. We solve this equation through 
iteration As our first approximation, we assume 

e(k + q) w e(k) which was used successfully by Mott and 



Jones to study metal conductivity for the reason that, 
when colliding with atoms, electrons tend to change mo- 
mentum rather than jump into orbits already occupied 
by other electrons pi. In addition to this, we are go- 
ing to see that e(k + q) ~ e(k) is actually the condition 
to have a superconducting gap. We consider a spheri- 
cal Fermi surface, so that both A and E depend only 
on e. Therefore our assumption e(k + q) ~ e(k) ensures 
A(k + q) sa A(k) and E(k + q) sa £l(k), which lead 
through Eq. M to 



E = Eq tanh 



E 



2k B T 



where 



En — 



E 

q 



(2) 



(3) 



which is just the value of E at T = 0: tanh(E/2k B T) -> 
1 when T ->• in Eq. (g). 

The iteration method was not used by BCS to solve 
Eq. ([l]). Instead, they assumed V = const, in that equa- 
tion, and converted the summation there into integra- 
tion, giving 
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N(0)V 



tanh 



E \ de 
2k B T J E 



(4) 



which is one of the key equations of the BCS theory, 
where A^(0) is the density of state near the Fermi sur- 
face, uj an unspecified phonon frequency. Further ap- 
proximation was applied to find the well-known ratio 
2A/k B T c — 3.5 J|]. However, it is interesting that E 
from Eq. (||) is also the exact solution of Eq. (|j), pro- 
vided that huiN(0)V = Eq: approximation is not nec- 
essary. This can be checked easily by simple inspection. 
Since tanh(£/2/c B T) < E/2k B T, we find from Eq. (§) 
T < En/2k B which means T c = Enj2h B and therefore 



2Eo_ 

knTr 



(5) 



which reduces to 2A/k B T c = 4 when e = 0. Therefore 
the ratio 3.5 must be replaced by 4, regardless of whether 
Eq. (Q) is solved via the original BCS approximation or 



1 



not. Indeed, the ratio 3.5 does not appear to have much 
advantage: we find from 23 superconducting elements in 
the literature [Q that on average 2Eo/kBT c — 3.8 with 
standard deviation 0.597. The deviation is 0.627 from 
the ratio 4.0, compared with 0.672 from ratio 3.5. 

Now we evaluate Eq in Eq. (||) quantitively. We have 



nk,q) = E 



2fcj,(q)A4?(k,q) 



[Mq)] 2 -Kk + q)-e(k)P 



(0) 



where oji is the phonon frequency, I identifies phonon 
branch (excluding transverse phonons, which do not in- 
teract with electrons in TV-processes) ||. Eq. (||) is 
slightly different from that in |l[ where phonon polar- 
ization is neglected. The matrix element 



TiN 


1/2 r 


.2Mw,(q)_ 


Jn 



M l (k,q)=q l ^773 / ^k +q> a(r)5V(r)V> k , CT (r)dr 

(7) 

measures the strength of electron-phonon interaction. 
Here ip is the electron wave function, a spin \ or j, M 
mass of an atom, N number of atoms in unit volume, r 
coordinates in real space, f2o a volume surrounding the 
atom, r its boundary, and SV(r) — V(r) — V(r ), V 
being the potential field. We define qi as the i-th com- 
ponent of [/q, U being the 3x3 unitary matrix found 
when solving the classical equation of motion for the 
atom. Mott and Jones found matrix elements when J7o 
is the Wigner-Seitz cell [§. We find Eq. (@) when V(r ) 
is constant (this defines flo in a natural manner). 

We use free electron energy to evaluate Eq. (^|) . This 
means that the Fermi surface is spherical, our assump- 
tion when solving Eq. (Q). As a result, the denomina- 
tor of Eq. (||) becomes 4e_F£ q [£ 2 — (C + cos6>) 2 ], e F = 
(h 2 /2m)\\i\ 2 is the Fermi energy (we study |k| near the 
Fermi surface), e q = (fr 2 /2m)|q| 2 , Sf = (m/2)vf /e_F, 
vi = cj;(q)/q| the sound velocity. In the Debye approx- 
imation Si = (Z/16) 1 / 3 6d/©_f ~ 10~ 3 in all supercon- 
ducting metals, where Z is the valency, Qd and Qf are 
the Debye and Fermi temperatures. It is apparent that 
V(k, q) > (condition to have an energy gap) holds 
in Eq. (§) onl Y whcn C + cosfl « 0, C = |q|/|2k|, 9 
being the angle between k and q, so that |k + q| 2 = 
|k| 2 + |q| 2 + 2|k||q|cos6> « |k| 2 . Thus e(k + q) « e(k): 
electrons change momentum but not energy in scatter- 
ing, the assumption by Mott and Jones ||. 

We substitute Eq. (^|) into Eq. (||) and replace the 
summation over q with an integration over (Air/tyk^ 
(volume of the first Brillouin zone), which exists in the 
sense of the Cauchy principal value (used by Kuper to 
verify the BCS theory) |^|, @], i.e. positive and negative 
contributions of V(k, q), if finite, are cancelled on a se- 
ries of spherical surface, the singular point ignored. We 
are entitled to do so, because Eq. (^|) is defined on a grid 
of k and q, which may not be in precise combinations to 
let V(k, q) = oo. We can also avoid such combinations 



by suppressing a few phonons with little physical conse- 
quence. This principal value varies little among phonon 
branches, allowing us to use J^i <1? = q*f^*t^q = |q| 2 {U 
unitary) to simplify the result. We use the expression for 
metal resistivity ||] to calibrate 6V in Eq. (0) and find 



E 



he 2 
kBT p 



■qripv 



(8) 



where e and n are electron charge and density, respec- 
tively, ks the Boltzmann constant, p the resistivity at 
temperature T p (not necessarily low), v = ksQo /hkc 
the Debye sound velocity, ko being the phonon cut-off 
wavenumber, and 
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(4Z)- 1/3 C 2 dC 



(42)- 



F 2 (x)( 3 d( » 1 
(9) 

Here F(x) — 3(a;cosx — sina;)/a: 3 is the overlap integral 
function, x = 3.84a 1/3 Z 1/3 C, a = Ntt /n the fraction 
of fio m a primitive cell, and f2 the unit volume. We 
assume |q|/|2k| = C < (4Z)- 1 / 3 < 1, because Eq. (§) 
arises from a canonical transformation where opera- 
tor commutation requires q ^ ±2k. 




FIG. 1. (A) Extended solution of the BCS self-consistent 
equation, E^ = E or — e, when — E < e < E or e < — E, 
respectively. (B) A = (E 2 — e 2 ) 1 / 2 which varies dramat- 
ically when — E < e < E. A « E only when e«0. 

It is apparent from Eqs. @j and (||) that E and Eq 
are constant. In contrast A = (E 2 — e 2 ) 1 / 2 varies dra- 
matically when — E < e < E. Beyond this range Eq. ([|) 
does not have non-trivial solution in first iteration. This 
justifies the approach of BCS to integrate Eq. (|l|) only in 
a thin layer across the Fermi surface Q . However Eq. (Q) 
has an extended solution (Fig. 1) 



£k = 



E 
— e 



-E < e < E 
e < -E 



(10) 



where E^ = — e arises from the trivial solution A = of 
Eq. (j|). The pair occupancy (Fig. 2) 



e 



(ii) 
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varies linearly when —E < e < E, but equals 1 when e < 
—E. This also justifies the BCS approach to integrate 
Eq. (Q) only in a thin layer: otherwise we have /ik > 1 
or < 0. The over-all probability of occupancy at T > 
is given by /ik(l — 2/k). Here 



/k = l 



1 



exp 



_Bk_ 



(12) 



is the probability of excitation which, unlike the Fermi- 
Dirac distribution, drops towards the interior of the Fermi 
sea (Fig. 2) where electrons are apparently more difficult 
to excite. Curves in Figs. 1 and 2 all have a kink at 
e = — E, which is likely to be smoothed out after further 
iterations. 




FIG. 2. (A) Electron pair occupancy, and (B) schematic 
of excitation probability. 

In the original BCS paper £7 k = (A 2 +e 2 ) x / 2 is known 
as the dispersion law In Eq. ( |l0| ) this dispersion 
law is reduced to the normal law once e < —E, consis- 
tent with the physics that no electron-phonon interac- 
tion takes place deep inside the Fermi sea. In the view 
of BCS the energy gap is measured by the minimum 
value of Ek Q . It is clear from Fig. 1 that E fits this 
definition. In comparison A appears to be a parame- 
ter to measure the electron-phonon interaction, which 
weakens dramatically with increasing distance from the 
Fermi surface, and cuts off when |e| > E. Furthermore, 
E in Eq. (||) matches experimental values of the su- 
perconducting gaps of metals reasonably well, despite 
that on surface the combination of various fundamen- 
tal constants in that equation appears to be a very long 
shot. Apparently, a good superconductor must have nu- 
merous free electrons (large n) scattered frequently by 
atoms (large p) moving quickly to facilitate pairing (large 
v). When a = 1, Eq. (|) yields 2E Q = 2.2, 18 and 
27 for Cd, Ta and Nb (1.5, 14 and 30.5 experimentally, 
in 10" 4 eV), which are of the right order, although over 
and under-estimations are possible. On average Eq. (|^) 
yields 2E = 15.7 for Zn, Cd, Hg, Al, Ga, Tl, Sn, Pb, V, 
Nb, Ta and Mo (11.3 experimentally). 

Although accepted by many, A = const, is from ap- 
proximation, employed to ease the mathematical diffi- 



culty to evaluate the BCS theory, and therefore not as 
fundamental as the theory itself. In fact A = const, leads 
to several conclusions that have to be justified. For ex- 
ample E — (A 2 + e 2 ) 1 / 2 increases when e > 0. It is not 
clear why pairs above the Fermi surface are more difficult 
to break, despite their higher energy. The dispersion law 
£ k = (A 2 + e 2 ) 1 / 2 reduces to the normal law only when 
|e| >> A. It is not clear why this normal law is distorted 
inside the Fermi sea, where no electron-phonon interac- 
tion is supposed to take place. Furthermore A = const. 
leads to Eq. (Q) which is not easy to be compared with 
experimental data. 

In conclusion the BCS theory has significant power of 
prediction when evaluated via the well-established method 
of iteration. This could be of current interest, because 
some novel theories for high T c superconductors resemble 
the BCS theory quite closely || ||] . These theories are of- 
ten phenomcnological, but nonetheless detailed enough 
for one to evaluate the matrix element Therefore a 
formula similar to Eq. (||) might arise, which could be 
useful when comparing the theory with experiments. 
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